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Abstract 

The inclusive cross-section for production of a jet with a given transverse momentum off a heavy 
nucleus is derived in the BFKL framework with a running coupling on the basis of the bootstrap 
relation. The cross-section depends on the same three different coupling constants as the total cross- 
section unlike the cross-section for gluon production derived in the dipole approach. 

1 Introduction 

Attempts to introduce the running coupling into the BFKL dynamics have a long history. As early as in 
1986 L.N.Lipatov introduced the running coupling by using an approximate form of the BFKL equation 
and a semi-classical approach [T]. In his derivation a single running coupling as{r) appears with r —> 0. 
As a result he found that the cut in the complex angular momentum J transforms into a sequence 
of poles, which condensate in a certain region depending on the total pomeron momentum. Later in 
our papers we introduced the running coupling into the BFKL equation by imposing the so-called 
bootstrap condition necessary for the fulfilment of unitarity mm- In this approach three different 
running couplings appear. The structure of the singularities in the J plane turned out dependent on 
the assumpsions about the uncontrolled low energy behaiour of the coupling. Remarkably it turned out 
that the bootstrap method correctly reproduced the running of the coupling in the reggeon interaction 
both in the forward and non-forward directions mm explicitly calculated in PIZIIH]. Still later, 
with the advent of the dipole picture and construction of the Balitski-Kovchegov equation for the 
pomeron in a heavy nucleus the running coupling was introduced by explicitly taking into account 
quark-antiquark loops in the evolution of the gluon density p uni ttn [n]. Again three different 
coupling constants appeared in the final equation. Remarkably this procedure turned out to be fully 
equivalent to the bootstrap approach, which leads to the same results with much less labour m- 
Finally a few years ago the dipole approach was generalized to the case of inclusive production off a 
heavy nucleus M- There at the leading order the number of different coupling constants proliferated 
up to seven. Still the authors made a conjecture for the form of the final inclusive cross-section, in 
whih the number of different coupling constants is reduced to three with one of them depending on 
the collinearity of the final jet components Acou, that is essentially on the experimental setup. 

In this paper we study the inclusive cross-section off a heavy nucleus within the bootstrap approach. 
It enables us to introduce the running coupling keeping the same number of them (three) as for the 
total cross-section. All the coupling constants remain depending only on the QCD scale A, with no 
dependence on the experimental conditions. Actually this approach allows to obtain only the inclusive 
cross-section for the production of a jet with a given transverse momentum, without specifying its 
other characteristics. 

Note that this observable is different from the one calculated in |14] . where the jet (or a ’’particle’) 
not only had a fixed transverse momentum but also had an upper bound on its collinearity Acou- The 


latter was assumed to be small and only the singlular terms terms at Kcoii —?> 0 were kept. In the 
case of hadron production it appears that Kcoii would become the factorization scale between the hard 
process and fragmentation function. In our approach seems to be effectively taken to infinity, 
making both the observable and its calculation very different. 

In any case our cross-section has a well determined physical meaning and can well be related to 
experimental observations. So we believe that it is worth studying, especially since it turns out much 
simpler than the conjectured cross-sections in m, which are more specific and so dependent on the 
experimental resolution. 

In our derivation we use the well-known two-dimensional transverse picture of BFKL and J.Bartels, 
(see e.g. m) We invoke the relative coefficients between different configurations in assumed corre¬ 
spondence with the standard AGK rules, which is confirmed by calculations with a fixed coupling. 


2 Triple-Pomeron vertex 


As mentioned in the Introduction, in this paper we follow the idea to introduce a running coupling 
via the bootstrap mm- Derivation of the triple-pomeron vertex in the limit Vc —>■ oo then goes as 
presented in [HllIT! for the fixed coupling case. The same derivation for the running coupling case 
is presented in | 13 ] . which we briefly recapitulate here. 

Basic formulas for the introduction of a running coupling via the bootstrap condition consist in 
expressing both the gluon trajectory uj and intergluon interaction in the vacuum channel V in terms 
of a single function rj(q) of the gluon momentum, which then can be chosen to conform to the high- 
momentum behaviour of the gluon density with a running coupling. Explicitly 
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In these definitions it is assumed that 91 -|- 92 = gi + g2 = g- Note that the unsymmetric form 
([2]) assumes that the initial pomeron with momenta q'l and 92 is ’’amputated”, that is multiplied by 
r]{qi)'i]{q2)■ This factor appears explicitly in the denominators in ([ 2 ]). For arbitrary 77(9) the following 
bootstrap relation is satisfied: 
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The fixed coupling corresponds to the choice 
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Then one hnds the standard expression for the trajectory 07(9) and 

where Vbfkl is the standard BFKL interaction. Note that the extra 27 r in the denominator corre¬ 
sponds to the BFKL weight l/( 27 r)^ in the momentum integration, which we prefer to take standardly 
as l/( 27 r)^. 

From the high-momentum behaviour of the gluon distribution with a running coupling one finds 
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As to the behaviour of rj{q) at small momenta, we shall assume 


r?(0) = 0, 


( 7 ) 


which guarantees that the gluon trajectory (jj{q) passes through zero at g = 0 in accordance with 
the gluon properties. The asymptotic Q and condition ([ 7 ]) are the only properties of T ]{ q ) which 
follow from the theoretical reasoning. A concrete form of r ]{ q ) interpolating between ([ 7 ]) and dS]) may 
be chosen differently. One hopes that the following physical results will not strongly depend on the 
choice. 

Our old derivation in m of the triple-pomeron vertex was actually based on the property ([Tj) 
obviously valid for (j 4 ]), the bootstrap relation and the relation between the Bartels transition vertex 
for a fixed coupling constant and intergluon BFKL interaction Vbfkl (Eq. ( 12 ) in [ 16 ]) 

1E/4"3(1,2,3|T,3') = VbfklMI' - 1,3') - Vbfkl{12,2,\1'^'). (8) 

Here and in the following we frequently denote gluon momenta just by numbers: 1 = gi, V = q'^ etc. 
Also we use 12 = qi + q2- All the rest conclusions were obtained from these three relations in a purely 
algebraic manner. 

If we define the transition vertex for a running coupling by a similar relation in terms of the new 
intergluon vertex V, Eq. ([ 2 ]), 

1E2^3(1, 2,3|T, 3') = E(2,3|1' - 1, 3') - E(12, 3|l'3'). (9) 


where V is given by ([ 2 ]) then all the results will remain valid also for arbitrary r]{q) satisfying ([ 7 ]) and 
thus for a running coupling, provided r]{q) is chosen appropriately. In this way in the momentum space 
one can find the the three-pomeron vertex with the running coupling m and write the corresponding 
Balitski- Kovchegov equation . It has a simpler form in the cordinate space. In the forward case at 
fixed impact parameter b the resulting Balitski-Kovchegov equation with the running coupling m is 
obtained as 
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where the two running coupling constants as{r) and Q;s(ri,r2) are defined as 
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and 

firi,r 2 } = J d^pfi{p)i{ri - p)^{r 2 - p). (12) 

where fj{r) is the Fourier transform of r]{q) and ^(r) is the Fourier transform of l/r]{q). Eq. (IlOp 
formally coincides with the running coupling equation obtained in the dipole aproach in pin]- 
However in our approach the couplings as(r’) and as{ri,r2) are not directly taken from the QCD. 
Rather it is function ^{q) that is taken from the QCD and the coupling constants are determined by 
it. 


3 Inclusive cross-sections 

Diagrammatically contributions to the forward amplitude for the collision of the projectile with the 
nucleus are identical to those which appear in the fixed coupling case. The difference is totally 
restricted to the new form of the reggeized gluon trajectory, interaction between the reggeons and the 
splitting vertex. Correspondingly the inclusive cross-sections will be obtained either by cutting the 
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Figure 1: Diagrams for the cross-section corresponding to emission from the upper pomeron (1-3) and 
from the triple pomeron vertex (4-6). Thick solid lines denote reggeons, thin lines denote real quarks, 
wavy lines denote real and virtual gluons. 


interaction in the uppermost pomeron (before all splitting) or by cutting the first splitting vertex. All 
the rest contributions will be eliminated by the AGK cancellations. 

However one should take into account that the inclusive cross-sections obtained in this manner 
do not refer to precisely gluon production. The form of the function r]{p) desribing the intermediate 
s-channel state with the single-loop /3 function includued implies that one has to take into account not 
only the single real gluon state but also states which contribute to this /S-function, namely the quark- 
antiquark and two-gluon states. So the inclusive cross-section obtained by fixing the intermediate s- 
state momentum p actually refers to all possible states having this momentum. It does not discriminate 
between contributions from gluons (single or in pairs) and (anti)quarks. So it is rather the inclusive 
cross-section for emission of a jet with a total transverse momentum p. The diagrams for this cross- 
section are illustrated in Fig. [TJ We think that this quantity has a well-defined physical meaning and 
is accessible for experimental observation and so worth theoretical investigation. In all the following 
we study precisiely this generalized inclusive cross-section for jet production. 

4 Jet production from the pomeron 

This part of the total cross-section is calculated in full analogy with the fixed coupling case. If only a 
single pomeron is exchanged between the projectile and target then we have the impulse approximation 
contribution 

IiY,y,p)=^-^^ = ^ r]{k)r]{p-k)P{Y-y,p-k)Piy,k). (13) 
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Here it is assumed that the ”semi-amputated” forward pomeron wave function 4>{y,p) = vip)P{y,p) 
in the momentum space satisfies the equation 

(M) 

With a fixed coupling P{y,p) obviously has order so that I has order g^, which takes into account 
two impact factors each of order g'^ 

For a nucleus target the pomeron coupled to the target has to be substituted by the solution of the 


























Balitski-Kovchegov equation (11011 at fixed impact parameter h (transformed to the momentum space): 
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5 Jet production from the vertex: generalities 

The three-pomeron vertex has a fixed rapidity and does not include evolution. This makes it possible 
to study the contribution to jet emission from the vertex in the lowest order of perturbation theory, 
that is for the target consisting of only two centers. It also allows to simplify treatment choosing for the 
projectile and target quarks, modeling the pomeron exchanges by colorless double reggeon exchanges. 
The inclusive cross-section will be obtained from the diagrams for the forward scattering off two 
centers, in which reggeon interactions and splittings are described by functions V and W, given by 
Eqs. (j 2 ]) and ([9]) with the running coupling. All the diagrams may be divided into two configurations 
depending on the way the four final reggeons are combined into pomerons. The diffractive one includes 
diagrams with two consecutive colourless exchanges (Figs. [2][4]). The non-diffractive configuration 
includes diagrams with parallel colourless exchanges (Figs. [MZl) when one of the colourless pair of 
reggeons is enclosed in the other. In these diagrams vertical lines denote reggeons, solid horizontal 
lines denote projectile and target quarks and wavy lines denote s-channel gluons. Diagrams in which 
the colourless pairs partially overlap do not give contribution to the inclusive cross-section and we do 
not show them. Note that the number of reggeons coupled to the projectile may vary from two to 
four. However, as shown long ago, all contributions reduce to that of a colourless pair of reggeons 
forming the initial pomeron (amputated with the forms ([2]) and ([9])). If one subtracts the contribution 
from the so-called reggeized term (see e.g. m) or alternatively from the Glauber initial condition in 
the dipole formalism) the the rest gives the contribution from the triple pomeron vertex, which is our 
goal 

To obtain the inclusive cross-section for the production of a jet with momentum p one has to 
fix function r](p) in intermediate states in the s channel. These intermediate states are obtained by 
cutting the diagrams in the s-channel. Different cuts may pass through one of the targets (single cuts, 
S), through both targets (double cuts, DC) or do not pass through targets at all (diffractive cuts, 
D). In the diffractive configuration only diffractive and single cuts are possible. In the non-diffractive 
configuration only single and double cuts are possible. According to the AGK rules the relative weights 
of the contributions from diffractive, single and double cuts are 1 ; — 1 ; 2 . 

We denote the reggeon momenta of the first final pomeron as qi and q 2 and of the second as q^ 
and ( 74 . In the forward direction q 2 = —qi and q^ = —q^. So the final pomerons carry momenta qi 
and q^. 


6 A. Diffractive configuration 

We divide all contributions according to the number of initial reggeons (R): A 2 , A3 and A4 for two 
three and four initial reggeons respectively. 

6.1 A4: 4R^ 4R 

We have 6 diagrams shown in Fig. [2j They contribute either in D or S configurations. Due to assumed 
relation between D and S diagrams 2 and 5 are ancelled. 1 and 6 give only S. 4 contains only D. Finally 
3 contains D-|-2S=-D. Due to the properies of the pomeron SDi = SDq = 0. So we are left only with 
3 and 4. Golour factors: 

^3 = ^4 = 

So taking (—H) for the interaction we find 


( 3 ) = -C5{-V{qi,q4\qi+p,q4-p)) 


( 16 ) 







Figure 2: 4R—)• 4R in the diffractive configuration 



Figure 3: 2R—?■ 4R in the diffractive configuration 

(4) = C4{-V{-qi,-q4\ - qi + p,-q^ + p))- (17) 

Changing signs of qi and q^ does not change V. So in the end (3)+(4)=0 and the diagrams give 
no contribution. 

6.2 A2: 2R^ 4R 

There is a single diagram 2R^4R, shown in Fig. [3l R contributes in both D and S configurations. 
Colour factor is C = N^. The D contribution is 

Md =-CW{qi,-qi-qi,qi\p,-p). (18) 

The two S contributions give 

Msd = +CW{qi,-qi - qi,q4\qi +P,-qi - p) +CW{qi,-qi - ^ 4 ,^ 4 ! - qA+P,qi -p)- (19) 


6.3 A3: 3R^ 4R 

We have 4 diagrams which contribute in D or SD configurations, shown in Fig. |4] Note that Ci = 
C 2 = - 1/2 and Ca = (74 = + 1 / 2 . 

So the S contribution from ( 1 ) cancels the left S contribution from (3) and the S conttibution from 
2 cancels the right contribution from (4). Moreover the D contribution from (3) cancels the right SD 
contribution from (3) and the D contribution from (4) canceles the left S contribution drom (4). As 
a result the only remaining contribution is D contribution from diagrams (1) and (2) 

43 = 43d = -CiW{qi,-qi,-q4\p,-qi - p) - CiW{-qi,-qi,q^\ - qi +p,-p), (20) 

where we have taken into account that D contributions from (1) and (2) are equal. 
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Figure 4: 3R—>'4R in the diffractive configuration 



Figure 5: 4R—)■ 4R in the non-diffractive configuration 

6.4 The total diffractive contribution 

Suppressing the overall colour coefficient we have found 

A = A 2 d + A 2 S + A 3 D = -W{qi, -qi - q^, q^lp, -p) 

+W{qi,-qi - q4,q4\qi +P,-qi - p) + W{qi,-qi - q4,q4\ - q4+P,q4-p) 

-^W{qi, -qi, -q^lp, -q^ - p) - ^W{-qi,-q4,q4\ - qi +p,-p)- 

7 B. Non-diffractive configuration 

7.1 B4: 4R^ 4R 

Six diagrams with transitions 4R—> 4R are shown in Fig. [5j Since the double cut (DC) contribuition 
enters with coefficient 2 diagram (3) is cancelled. Diagrams 2 and 5 give the same contribution. The S 
in both cancels one half of the DC contribution, so that the total coefficient is reduced to 1. Diagram 
(4) is zero. Diagram 1 and 6 contribute only to the S contribution. So in the end the DC and S 
contributions are 

R4z5c = (2) + (5), BAs = -{!)-{6). (21) 

Colour factors are 

C2 = C5 = -\n^, Ci = c, = ^Nt 
So we get (again suppressing N^) 

BAdc = ^{y{-q4,-qi \ - q4+P,-qi -p), +V{qi,q4\qi +P,q4-p)), 


( 22 ) 
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Figure 6: 2R—>■ 4R in the non-diffractive configuration 
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Figure 7: 3R—)■ 4R in the non-diffractive configuration 


B4:s = ^(v{-q4,qi\ - q4+P,qi - p) + V{-qi,q4\ - qi +P,q4-p))- 


(23) 


7.2 B2: 2R-> 4R 

There is a single diagram for transitions 2—4 shown in Fig. [6l R contains DC and S contributions. 
The colour factor is C = N'^, so suppressing it 

B2dc = -2bF(-^ 4 ,0, -q^+P, -Qi + Qa - p) (24) 

and 

B2s = 2W{-q4,0,q4\ -q4+P,q4-p)- (25) 

In fact the S contribution gives zero: it does not depend on qi. 


7.3 B3: 3R^ 4R 

Four diagrams with transitions 3—T'd are shown in Fig. [7l We observe that the right S contribution in 
(3) cancels one half of the DC contribution in (3) and the left S contribution in (4) cancels one half 
of the DC contribution in (3). Thus we get 


B3dc = 2(1)_dc + 2(2)£)C' + (3)dc + {‘^)dc 


and 


B2,s = -(l)s - (2)5 - (3)s/ - (4)sr, 


where SI and Sr refer to left and right single cuts. 
The colour factors are 


Ci = C2 = \Nt, C^ = C4 = -\Nt 


So suppressing 


(26) 

( 27 ) 


B‘iDC = - W(-q4, qi, -qi 
+ + ^w(-q4,qi,q4lqi 


-q4 + qi +P,-qi -P) -h -W(qi,-qi,q4lqi +P,q4-qi - p) 
q4+P,q4-p) + ^W{-q4,-qi,q4\ - 94 +P, 94 - 9i - p) 




































and 


B 3 s = -W{-qi,qi,-qi\ - q^ + p, -p)-W{qi, -qi, q^lp, q^ - p) 
-^W{-q4,qi,qA\ - qA+P,qi + Qa- p) - ^W{-q4,-qi,qA\ - qi - qA+P,qA-p)- 

7.4 The total non-diffractive contribution 


B = BAoc + BAsdc + B2dc + B2s + B3 dc + B3s = 

^(v{-qA, -qi\ - qA+p,-qi -p) + y{qi,qA\qi +p,qA-p) 

+ {y{-qi^qi \ - qA+P,qi -p) + V{-qi,q4\ - qi +p,qA-p)) 

-2W{-qi,Q,q4\qi -qA+P, -qi +qA-p) 

-W{-qA,qi,-qi \ -qA + qi +P,-qi - p) - W{qi,-qi,qi\qi +p,qA- qi - p) 

+^{-qA,qi,qA\qi - qA+p,qA-p) + ^{-qA, -qi,qA\ -qA+p,qA- qi -p) 

+ ]^W{-qA,qi,-qi\ - qA+P,-p) + ]^W{qi, -qi,qA\p,qA - p) 

-ilT(-g4,gi,g4| - qA + p,qi+ qA- p) - ^{-qA,-qi,qA\ -qi- qA+p,qA-p)- 

8 Jet production from the vertex: the final form 

The found contributions to inclusive cross-sections from the two configurations are numerous and 
contain many terms. To analyze them it is convenient to introduce function 

F{q\p) = {q\p) = -P) = (-'?! - P)- (28) 

In terms of this function 


V{qi,q2\ki,k2) = {ki2\ki) - (gi|A;i) - (^21^2), 


(29) 


W{qi,q2,q‘i\ki,k2) = {ki2\ki) - (^121^1) - (^231^2) + (g' 2 |^i - qi), (30) 

where ki 2 = ki + k 2 etc. The analysis of the found inclusive cross-section in terms of function F{q\p) 
is performed in Appendix. It is important to take into account that all contributions which do not 
depend on the momenta of one of the final pomerons, qi or 54 , vanish, since the Pomeron vanishes 
when its two reggeized gluons are at the same spatial point. 

After some cancellations one finally finds the contribution 


A + B 
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{qiA\p) + {qiA\ -p) 


(31) 


With the hxed coupling constant one finds instead 

{A + B)B^ = 

TT 

It was shown that in the hxed coupling case half of this expression comes from the reggeized piece (or 
the Glauber initial condition in the dipole approach) |19] . By the same reasoning the contribution 
from the triple pomeron vertex is one half of (|31l) . As we see, to hnd our inclusive cross-section for 
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jet production with the running couling all the necessary change is the substitution of all momenta 
according to the rule 

^ (33) 

Using the known form of the inclusive cross-section in the fixed coupling case ([IS], see also [T9| for 
closer notations) one can immediately write the final expression for it with the running coupling as 


h{Y,y,p) = J - p)P{Y - y,qii - p)r]{qi4)<^>{y,qi)<^{y,q4). 

The total inclusive cross-section is the sum of (fT5]) and ([3^ : 

{2'n)^da 2N^ 


(34) 


I{Y,y,p,b) = 


(PqicPqi 


dyd?p(Pb p{p) J (27r)^ 
p{q(U - p)P{Y - y, 914 - p) \2r]{qi)<^>{y, qi,b)){2TTf5{qi) - p{qi4)<^{y, 9i)4>(y, 94 ) 


(35) 


9 Discussion 

Our final expression for the inclusive cross-section Eq. (I35p has a strong similarity with the one 
conjectured in [Hj). It contains three factors p, which can be put in correspondence with the three 
coupling constants depending on different arguments in m- However in our formula the arguments 
of functions p directly depend on the three momenta involved: that of the observed jet and two of the 
gluon distributions involved. In contrast, in the conjecture of m the argument of one of the coupling 
constants depends only on the assumed value of collinearity of the observed jet and the arguments of 
two others are complex and depend on all three momenta in a very complicated manner. However, as 
mentioned in the Introduction, the literal comparison of the two cross-sections is not possible, since 
in fact they refer to different processes: to jet production with a given momentum in our case and 
with additional restriction on jet collinearity in [Tl|. Besides the cross-section in |14| is after all 
only conjecture, whereas ours is more or less consitently derived from the bootstrap condition which 
demonstrated its validity for the total cross-section. 

It remains to be seen by practical applications to what extent this difference is felt in the actual 
inclusive cross-sections. To do this consistently one has to previously solve our equation for the 
unintegrated gluon density ([TO]) . With all its similarity to the currently used equation in the dipole 
picture, the actual values of the three running couplings involved are not identical, so that the already 
found solutions in the dipole picture cannot be directly used for our inclusive cross-section. We 
postpone this problem for future studies. 

In conclusion we stress that that our equations take in account terms of the order {a{q))'^ with a{q) 
taken in the leading order. Subleading terms of the relative order l/ln( 9 ^/A^) remain undetermined, 
since they correspond to the next-to-leading order in the running coupling. 


10 Appendix.Analysis of the found inclusive cross-sections 

In terms of F we find 


A = -(0 - (-94|p) - (-gil -p) + (-gi - q4lp-gi)) 

+ (o - (-q4lqi +p) - (- 91 I -qi-p) + (-qi - q4\p)) 

+ (o - (- 94 I - q4+p) - {-qi\q4-p) + {-qi - q4\ - p)) 

-q4\p) - {-qi - q4\ -q4-p) + {-qi\p - qi))^[{-qi\ -p)- {-qi - q4\ -qi+p) + {-q4\p)) 


1 r\ 









We take into account that terms which do not depend on qi or do not depend on §4 are integrated 
ober qi and q^ to give zero. Dropping these terms 


^ = -{-qi - 94b - 9 i) - (- 94 ki +p) + {-qi - q4\p) 

-{-qi\q4-p) + {-qi - q4\ -p) 

-^(-^1 - 94 I -q4-p)- \(-qi - qi\ -qi+p) 

Now non-diffractive contributions 

B = ^((-^4 - qi\ -q4+p)- {-q4\ - q4 + p) - (-gi| - qi - p) 

+{qi + ^4bi +p) - {qi\qi +p) - {q4\q4-p) 

+ {qi - q4\ -q4+p)- {-q4\ -q4+p)- {qi\qi - p) 

+{q4 - ^i| - qi +p){-qi\ - qi +p) - {q4\q4-p)) 

-2( - {-q4\qi -q4+p)- {q4\ - qi + q4 - p)) 

-((- 94 I -qi-p)- {qi - q4\qi - q4 + p) + {qi\qi + p)) 

-{{q4\qi +p)- {q4 - qi\q4 - qi - p) + {-qi\p)) 

+^{iqi\q4-p) - {qi - q4\qi - q4+p) - {qi + q4\q4 - p) + {qi\qi +p)) 

+ ^((- 91 1 -q4+p)- {-qi - q4\ -q4+p)- {q4 - qi\q4 - qi - p) + (-^ib)) 

+ ^((-^4| -p)- {qi - q4\ - q4+p) + bi| - p)) + ^[{q4\p) - {q4 - qi\q4 - p) + (-gi| - p)) 
-^((9i| -q4 + p)- {qi - q4\ - q4+p) - {qi + q4\qi + q4 - p){qi\p)) 

-^{i-qi\q4 - p) - i-qi - q4\ -qi- q4+p) - {q4 - qi\q4 - p) + (-gi| - p)) 

Again we remove terms which do not depend on one of the variables qi or (74. Dropping these 
terms 

B = |((-94 - qi\ - q4+p) + {qi +94^1 + p) + {qi - q4\ - q4+p) + {q4 - qi\ - qi + p)) 

-2( - {-qi\qi -q4+p)- {q4\ - qi + q4 - p)) 

-((-94I -qi-p)- {qi - q 4 \qi - 94 +p)) - ( - (94 - 9 i |94 - 9 i - p) + (94I91 +p)) 
+^(( 9 i |94 -p) - (91 - 94I91 - 94 +p) - (91 +94I94 -p)) 

+ |((- 9 i| - 94 +P) - (-91 - 94I - 94 +P) - (94 - 91I94 - 91 -P)) 

+ ^( - (91 - 94I - 94 + P)) + ^( - (94 - 91I94 -p)) 

-^(( 9 i| -94 +p) - (91 -94I -94+p) - (91 +94I91 +94 -p)) 

-|((- 9 i |94 -p) - (-91 - 94I - 91 - 94 +P) - (94 - 91I94 -p)) 



Next we simplify our expressions using the possibility to change 91^4 -H- —( 71,4 and g' 1^4 -f-)- q^^i. We 
reduce all terms to two basic structures (gi 4 |o) and (gi|a). Then we get 


^ = -{qi^qi +p)- {qi\qi+p) + [qi^p] - {qi\qi-p) + [qul -p) 
-\{qiAqi +p)- ^(914^1 +p) 

and 

B = {qiAqi +p) - 2( - (gi|g4 -p) - (qilqi+p)) 
-{qi\q4-p) + (^mI -p) + (quip) - {qi\q4+p) 
+^{{qi\q4-p) - {qi4\ -p) - {qi4\qi +p)) 

+^{{qi\q4 + p) - {qi4\qi +p)- {qi4\p)) - {qulqi +p) 
-\{{qi\q4+p) - {qi4\qi +p)- {qi4\p)) 

-^{{qi\q4-p) - iqi4\ -p)- {qi4\qi +p)) 

Note that 

( 914^1 -p) = {qi4\q4-p) = {qi4\qi+p) 

Let us find the coefficients for our basic structures 

(914 1 gi +p) 

A2:-l, A3:+l, B4:+1,B2: 0,B3:-l/4-l/4-l+l/4+l/4=-l, total=0. 

{qi\q4-p) 

A2:-l, A3:0, B4:0, B2:2, B3:-l+l/4-l/4=-l, total=0. 

{qi\q4 + p) 

A2:-l, A3:0, B4:0, B2:+2, B3:-l+l/4-l/4=-l, total=0. 

(^14 Ip) 

A2:+l, A3:0, B4:0, B2:0, B3:+l-l/4+l/4=+l, total =2. 

(914 1 -p) 

A2:+l, A3:0, B4:0, B2:0, B3:+l-l/4+l/4=l, total =2. 

This brings us to the final result (ISTI) . 
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